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Hybrid architectures integrating mesoscopic electronic conductors with resonant microwave cavities have a 
great potential for investigating unexplored regimes of electron-photon coupling. In this context, producing 
nonclassical squeezed light is a key step towards quantum communication with scalable solid-state devices. 
Here we show that parametric driving of the electronic conductor induces a squeezed steady state in the cavity. 
We find that squeezing properties of the cavity are essentially determined by the electronic noise correlators of 
the quantum conductor. In the case of a tunnel junction, we predict that squeezing is optimized by applying a 
time-periodic series of quantized 5—peaks in the bias voltage. For an asymmetric quantum dot, we show that a 
sharp Leviton pulse is able to achieve perfect cavity squeezing. 

PACS numbers: 73.23.-b,72.70.+m,42.50.Ct,84.40.Az 


I. INTRODUCTION 


Squeezed states of light IT] exhibit reduced noise below the vacuum level in one of their quadrature and amplification in 
the other quadrature. Their realization is a key step in the development of quantum communication. They are important tools 
for continuous variable quantum information protocols EEl where they serve as building blocks for generating non-classical 
states. Their enhanced sensitivity can also be used for quantum non-demolition measurements of position and force H. Easily 
produced in optical systems, squeezed states have been observed more recently in circuit quantum electrodynamics at microwave 
frequencies a, either as single-mode EE), two-mode squeezing Ena or as Einstein-Podolsky-Rosen states Qo). 

The parametric driving used so far in experiments is limited to a half-squeezed quadrature for a cavity mode because of the 
inevitable coupling to the external vacuum fluctuations m. This limit however does not apply to dissipative squeezing, in which 
one steers the environment to stabilize the cavity into a non-classical state. In this case, perfect squeezing can be achieved, at 
least in principle, with minimum uncertainty ifT^ [TTll . 

A recent development in the held of superconducting quantum circuits is the realization of hybrid systems in which a quantum 
conductor is coupled to a microwave resonator. These systems offer an appealing platform for investigating fundamental matter- 
light interactions with an experimental control on both the electronic and photonic parts HlHIl]- Experiments have been realized 
with metallic tunnel junctions connected to a resonating line mi , or with quantum dots, realized in carbon nanotubes, nanowires 
or two-dimensional electron gases EOl - Ehl embedded in high-finesse coplanar cavities. The interplay of electron transport and 
emission of photons can lead to an electronic-induced lasing state in the cavity lEzHlo], and more generally produce bunched or 
antibunched photons EiKia, and nonclassicality in the light emitted by a quantum conductor ESI El. Squeezed light emitted 
by a tunnel junction was recently demonstrated experimentally in the absence of a cavity Il38ti40l . 

In this paper, we describe dissipative squeezing of a cavity mode coupled to an ac driven electronic reservoir. The system, 
depicted in Eig. [T] is a quantum conductor coupled to a microwave resonator. In addition, a classical bias voltage, with an 
oscillating part at twice the resonator frequency, is applied to the conductor. The conductor plays the role of a nonlinear 
environment: photons from the ac modulation are broken into pairs and transmitted to the cavity, thereby producing squeezing. 
We show that the amount of cavity squeezing is determined solely by current noise fluctuations in the conductor. We focus our 
attention on ac excitations which optimizes squeezing. In the case of a tunnel junction, we And that the best solution consists 
of periodic and quantized voltage peaks occurring in phase with the compressed quadrature. We also discuss how squeezing 
improves with the number of harmonics in the ac signal. Eor a quantum dot, we identify the conditions of optimum squeezing: 
asymmetric coupling to the leads, a narrow single-level resonance, a far-detuned single-level energy and a dc bias voltage 
matching the resonator frequency. In addition, we show that a Leviton pulse results in a vacuum squeezed state with minimal 
uncertainty. Perfect squeezing is approached by narrowing the width of the voltage pulses in the Leviton, again in phase with 
the compressed quadrature. 
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FIG. 1. (a) Schematic of a resonant cavity realized in a superconducting transmission line. The center line is coupled to external electromag¬ 
netic modes (ports) via capacitive gaps delimiting the cavity. The output field of the cavity can be measured in these ports. A dc voltage can 
be applied to the center line at a voltage node in order to preserve the high-Q factor of the cavity EHED- (b) Focus on the galvanic coupling 
between the cavity and a quantum conductor. The electric field (arrows) is spatially constant at the scale of the conductor. The conductor 
is lead-contacted on the upper side to the center line and on the lower side to the ground plane, (c) The quantum conductor can be a tunnel 
junction or a quantum dot. The upper lead is connected to the center line where the voltage potential is the sum of dc [V] and ac [14c(t)] parts, 
as well as the cavity quantum field oc a -f a^. 


II. TUNNEL JUNCTION 

We start by considering the case of a tunnel junction for the quantum conductor. The voltage bias across the junction is the 
sum of a driven classical part applied on the upper lead, and a quantum part associated to the cavity field. Using a convenient 
choice of electromagnetic gauge, it simply dresses the electron tunneling operators 

iTT(f) (1) 

where the operator T transfers one electron from the upper to the lower lead while does the opposite ll4^|43l . is the total 
semi-classical phase accumulated during a tunneling event, decomposed as if{t) = {eV/h)t + (j)g,c{t)+i{g/u!o){a^ — a), with the 
ac phase (padt) = (e//i) /* dt' Vac{t'). The last term is the quantum part and g measures the junction-cavity coupling strength. 
The cavity Hamiltonian is reduced for simplicity to a single mode, iTcav = huj^cdd with the cavity annihilation operator d. 
Eq. Q contains both the excitation of the cavity state by electron tunneling events and photo-assisted transport phenomena 
triggered by the ac modulation iniia. H'Yd) is thus responsible for an exchange of energy between three sub-systems: the 
cavity, the ac classical field and lead (free) electrons. 

We assume weak junction-cavity coupling and therefore neglect the backaction-induced change in electron tunneling resulting 
from the cavity. We thus set p = 0 to examine the current fluctuations of the tunnel junction and latter reinstate a finite g when 
considering the dynamics of the cavity field d. In the presence of ac voltage modulation, the photo-assisted noise properties of 
the tunnel junction are characterized by the correlator (at p = 0) ll46ll 

+00 

S'n(wi) 271(5 (wi + 072 — 2nwo), (2) 

n— — oo 


where I (w) is the Fourier transform of the Heisenberg current operator I (t) of the junction. The n = 0 term gives the stationary 
part of the noise, meaning that, when Fourier transformed with times ti and t 2 , it depends only on the time difference U — t 2 
and not on the mean time i = (ti + f2)/2. So(uj) is called the absorption (emission) noise of the tunnel junction for w > 0 
(a; < 0) iZl. It governs the rate of energy transfer between the junction and its environment, here the cavity, via single photons 
of energy h\uj\. 

As shown in the Supplementary Note 1, the different noise terms can be calculated to leading order in the tunneling strength 
2S'„(u;) = ^ [cn'cd+nS{to + eV/h + 2n'wo) + - eVjh - 2n'u;o)] (3) 

n' 

in terms of the (unsymmetrized) equilibrium Johnson-Nyquist noise of the tunnel junction S{io) = (2w/i?T)(l— 
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where Rt is the junction dc resistance. The coefficients Cm are defined by the Fourier expansion 


c(t) _ 


= E' 

m^Z 


2imu}ot i 
7,C C 


(4) 


of the ac phase c/iac- is the overall phase of the ac signal. For the lead connected to the voltage Vac{t) = {h/e)(pac{t), Cm gives 
the probability amplitude for an electron to absorb m energy quanta from the classical ac field when m > 0. m < 0 describes 
correspondingly photon emission to the ac field ll48l . In the particular case of a sinusoidal excitation, 14c(f) = Vi cos{2ujot + (p), 
these coefficients are written in terms of Bessel functions 

The energy of absorbed photons, 2mhujQ, can be either used in exciting energetic electron-hole pairs in the conductor, or 
transferred to the cavity. The non-stationary noise terms Sn^o in Eq. (|^ oscillate with the mean time t. They do not conserve 
energy and can provide n quanta of energy 2hu!Q (or absorb if n < 0) to the cavity. We will return below to the physical 
significance of these terms when analysing the cavity stationary state and squeezing effects. 

Now that we have detailed the possible transfers of energy between the ac-excited tunnel junction and its environment (the 
cavity), we study the cavity evolution under the dissipative influence of the electrons. Assuming weak junction-cavity coupling, 
we expand Ht to first order in g/ujQ, Hr = -f iXI{a^ — a) with the coupling constant A = tig/{euJo). The cavity evolution 
is described by a Heisenberg-Langevin equation 


d + iujQd + —a = XI{t), 


( 6 ) 


justified, either by an input-output calculation ||4^ detailed in the Supplementary Note 2, or by a Keldysh path integral formu¬ 
lation, discussed in appendix]^ assessing the Gaussian character of current and cavity field fluctuations. The tunnel junction 
current I{t) in Eq. (|^ plays the role of a quantum noise term, with fluctuations characterized by the correlator of Eq. (|^. In de¬ 
riving this equation, we have neglected the intrinsic (bare) damping of the cavity kq, assuming that the cavity dissipation caused 
by the electrons dominates. The corresponding damping rate k = A^[S'o(a;o) — 5'o(— wq)] balances absorption and emission 
noises, since absorption (emission) noise corresponds to photon loss (gain) from the cavity. Our calculation is also based on the 
use of the rotating-wave approximation where rapidly oscillating terms are averaged to zero. The validity of this approximation 
is controlled by the smallness of k/ujq and is consistent with our first order conductor-cavity coupling and with the absence of 
cavity backaction. 

The first order differential Eq. can be solved straightforwardly in time or frequency space, and yields the steady-state 
correlation functions for the cavity field d. We find for the anomalous correlator, using that k <C ujq, 


X^Si(ojq) _ S'ljwo) 

K S'ojwo) - 5'o(-a;o)' 


(7) 


This result can be given a physical interpretation: (wq) describes the coherent emission by the junction of a quantum of energy 
2Hujq to the cavity. This energy quantum breaks into a pair of cavity photons thereby contributing to the (a^) correlator. This 
effect is limited in the denominator by the rate at which cavity photons are absorbed by the tunnel junction. In the same way, the 
number of cavity photons (a^a) = X‘^So{—u!o)/k is unsurprisingly governed by the electronic emission noise. 

We now investigate field squeezing more precisely and introduce the two quadratures Xi = — de“^^‘^) and 

X 2 = where we use the same phase tp as in Eq. Q. Their variance is readily obtained 

, ^2 ^ En \CnTCn+l\^Sei,eV/n+ {2n+l)uJo) 

En(|cnP-|c„+iP)5o(ey/S+(2n + l)cuo)’ 


where we introduced the even Sg and odd So parts of the Johnson-Nyquist noise S. The two quadrature fluctuations are thus 
sensitive to the electronic temperature T, the dc bias voltage V and the pulse shape of the ac signal. The squeezing mechanism 
is optimized by taking the limit of vanishing temperature and by setting el^ = Huiq. 

We first consider a single-tone driving of the tunnel junction, I4c(f) = Vi cos(2wot + ^p), the experimentally most acces¬ 
sible situation. The photo-assisted coefficients c„ are then given by Bessel functions, as detailed in appendix]^ A numerical 
minimization of AXf in Eq. ([^ gives AXf = 0.618 for eVi = 0.706 x 2hoJo, with AA| = 1.864. This optimal squeezing 
value coincides exactly with the squeezing of the emitted light predicted and measured in Refs. see also the more recent 

Ref |2Q1. This is explained by noting that the zero-temperature cavity damping k/A^ = 2f>wo/RT (the denominator in Eq. <0) 
is constant for a tunnel junction, regardless of the bias voltage shape. This independence no longer holds at finite temperature or 
in the case of a conductor with a non-linear I-V characteristic. 
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Number of Harmonics 


FIG. 2. Minimum variance of the squeezed quadrature (left) and Heisenberg uncertainty product (right) as a function of the number of 
harmonics in the ac signal applied to a tunnel junction. The variance AXi of Eq. is minimized numerically for the ac signal Vac (t) = 
Vj cos{2jujot) where n is the number of harmonics. Optimal squeezing is obtained with T = 0, eV = hujQ. For comparison, the 
variance has been computed for a Leviton pulse with the result = 1 — 2a + 2a^. The optimized Leviton with a = 1/2 is shown in the 
last column. 


We turn to an ac modulation with the same fundamental frequency 2a;o but including higher harmonics El. Fig.|^shows the 
improvement in squeezing by adding more and more harmonics while AX| is further amplified. Considering a general 
periodic signal, we find analytically, as shown in appendix]^ that the minimum value AA^ = T/tt^ = 0.405 is reached when 
Cn = (l/7r)(n + l/2)“^ for n S Z, in agreement with a numerical minimization. The corresponding ac phase across the junction 
is aperiodic piecewise linear function 0ac(f) = 7r/2 —wof forf g]0, 7r/a;o[, with a jump discontinuity of tt at f = 0 and multiples 
of tt/wq. Adding the dc voltage V, we find for the optimal voltage applied to the junction a series of 5-peaks, 

iGZ ^ 

It is useful to give an intuitive classical picture for this squeezing optimization; the bias potential 14pt(f) acts on the conductor 
specifically at times where the amplitude of the squeezed quadrature is maximum and the other quadrature vanishes. This is in 
fact a strong perturbation, the emission and absorption noises are infinite and the second variance AA| = -l-c». 

Despite its linear current response, the tunnel junction is able to squeeze the cavity state down to 40% of the zero-point level. 
This is because squeezing is not governed by the current itself but by current fluctuations, and the noise of a tunnel junction is 
only a piecewise linear function enabling rectification lf38l . Eq. ([^ nevertheless suggests that better squeezing can be achieved 
by using a genuine non-linear system. 


III. ASYMMETRIC QUANTUM DOT 

We consider a quantum dot for the conductor embedded in the cavity. The situation where the dot is symmetrically coupled 
to the two leads, discussed in the Supplementary Note 3, is not optimal for squeezing. We thus focus on the asymmetric case 
where the upper lead is more weakly coupled to the quantum dot than the lower lead. In this case, the voltage drop from the 
central strip to the ground mainly takes place at the upper dot-lead tunnel contact. The first order cavity-conductor coupling is 
then of the form iXIu{a^ — a) where Ijj denotes the electrical current of electrons incoming from the upper lead. The coupling 
to the lower lead current is neglected. 

In practice, for quantum dot geometries, it may be important to also take into account the coupling of electronic transport to 
phonons. Following Ref. Il52l . it would lead to tunneling processes involving the excitation of phonon-photon pairs, degrading 
the quality of squeezing. Such study is nevertheless beyond the scope of this work and we neglect electron-phonon coupling in 
what follows. 

The analysis developed above for the tunnel junction can be essentially carried over to the quantum dot, with lu replacing I 
in the Heisenberg-Langevin Eq. The noise properties of the quantum dot are derived using scattering theory as discussed in 
appendix [D] We retrieve noise factors similar to Eq. Q, 

2Sn{uj) = ^ [c„/c*,+„.S+(a; -F eV/h -F 2n'wo) + - eVjh- 2n'wo)] 

n' 


( 10 ) 
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FIG. 3. Squeezed variance for an asymmetric quantum dot driven by a single-tone excitation I4c(f) = Vi cos(2a;of + v) versus the dot 
single-level energy, as given by Eq. ([§]l and Eq. (Uj. AXj is computed at T = 0, el^ = fituo and V = O.Icuq. For each £d, Vi is tuned to 
minimize AXi. The first local minimum AAi = 0.358 is reached for = 1.54&uo and eV\/{2fujJo) = 1.16. Asymptotically at large £d, 
the minima follow the curve 0.436(£d/2fta;o)~^^^ (in dotted line) and the optimal ac amplitude follows eVi = £d- 


here involving two different equilibrium noise terms 


S±{u:) 


2e'^TuT 


J{e-Huj) [!-/(£)] 
{eTedf + {hr/ 2 y^ 


( 11 ) 


with the Fermi function/(e) = The broadening F of the dot single energy level, denoted can be decomposed 

according to its coupling to the upper and lower leads F = Fjy -F F^ with Fj/ <C F^. The Lorentzian form in the integrand 
of Eq. ( [TT| l describes the Breit-Wigner resonance for transmitting electrons through the dot jT). S+{uj) describes electron-hole 
excitations with energy hu, where the electron, with energy e, tunnels from the lower to the upper lead and has to meet the 
resonance condition of the dot single-level. S- is the same but with hole tunneling. Proceeding with the calculation of cavity 
properties based on the Heisenberg-Langevin equation of motion, we retrieve the quadrature variances of Eq. ([^ if we define 
the even/odd parts as 2S^/o{u}) = S+{u}) ± S-{—oj). 

We have studied numerically the minimization of the variance AXf. Quite generally, squeezing optimization requires zero 
temperature and the four-wave mixing condition eV — huo l38l, as well as a sharp resonance, F ujq. In this regime, for 
Ed > 0, S- becomes negligible and S+ in Eq. is either constant for n > n-th = Zd/ (2lA;o) — 1, or vanishingly small below 
this threshold. Se and So simplify in Eq. ([^ and the summation involves only values of n above the threshold nth- 

In the case of a single-frequency ac modulation, Eig. shows that the squeezed variance displays a series of local 

minima, where the values of the minima decrease with Ed- Large single-level energy Ed however implies stronger power in the 
ac excitation signal in order to meet the Breit-Wigner resonance condition. In practice, this requires an ac signal amplitude Vi 
close to Ed such that electrons can tunnel through the dot. Perfect squeezing AXf = 0 is reached at very large Ed but only with 
a weak power law of coefficient —1/3. 

Alternatively, a vacuum squeezed state can be reached in the cavity by using a Leviton ac signal with the fundamental fre¬ 
quency 2uJo. Leviton pulses were originally 1541 proposed as voltage excitations designed to transfer a finite number of electrons 
through a coherent conductor with minimal noise, in analogy with coherent states minimizing quantum-mechanical uncertainty. 
They consist of sums of Lorentzian pulses with unit flux each. A unit flux represents the attempt to transmit a single electron. 
Mathematical details are briefly reviewed in appendix[A| A Leviton can be periodized by having an infinite train of evenly spaced 
Lorentzian pulses l48ll55l . Leviton pulses have recently been synthesized and used to perform Hong-Ou-Mandel electronic ex¬ 
periments 15^ and electron quantum tomography El. 

The use of a Leviton pulse for squeezing is natural. A Leviton train with periodic short pulses addresses specifically one 
quadrature (the one to be squeezed), while producing a minimal disturbance (noise) on the quantum conductor. Taking the limit 
of a very sharp resonance F /ojq —> 0 and the bias voltage eV = hwo, one obtains for the squeezed variance 


AXf 


EZo(cl-cl^^) 


I — a 
1 -F a ’ 


( 12 ) 


for 0 < Ed < 2hujo, and a minimal Heisenberg uncertainty AXi AX 2 = 1 reflecting the minimized perturbation by the Leviton 
compared to other types of ac excitation, a = is a parameter related to the width r of each Lorentzian pulse. A Leviton 

pulse is thus able to produce in optimal conditions an ideal squeezed state with arbitrary compression. 
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FIG. 4. Squeezed variance for the asymmetric quantum dot as a function of the Leviton pulse excitation parameter a. The continuous line 
was computed using Eq. and Eq. (Uj at zero temperature, eV = Twjq, E/wo = 0.1, Ed = Oficuo and the Leviton Cn coefficients (see 
appendix]^. It is compared to the asymptotical expression (1 — a)/(l + a) (dotted line) valid for an infinitely sharp resonance F —>■ 0. Note 
that the limit a = 1 corresponds to a dc voltage eV^ = —hujo producing no squeezing AXj = 1, indicating that the limits E/cuo —>■ 0 and 
a —> 1 do not commute. 


As shown in Fig. squeezing saturates at finite damping E. For a reasonable damping rate F/wq = 0.1, we find that the 
squeezed variance can still be reduced down to AXf = 0.075 for a = 0.91. A Leviton excitation with parameter a close to 1 
exhibits sharp voltage pulses of width r ~ (1 — a) /{2ujq). Similarly to the tunnel junction, it corresponds to concentrate short- 
time pulses of high voltage when the amplitude of the squeezed quadrature is maximum while the other quadrature vanishes. 


IV. CAVITY READOUT 

Cavity squeezing could be evidenced by an in-situ measurement using a qubit and its anisotropic radiative decay ll58l . It is 
also possible to demonstrate squeezing in the cavity by measuring the output field. So far, we have neglected in our discussion 
the coupling of the cavity to the external electromagnetic modes, see Fig.[T] They introduce an additional damping rate kq- As 
detailed in the Supplementary Note 2, this external damping adds two terms in the Heisenberg-Langevin Eq. (|^, which now 
reads 

a + fcuod + = -v^^in(f) + A/(f). (13) 

bin{t) is the input field, it describes the quantum state of incident photons on the cavity. In the absence of an input drive, it 
corresponds to vacuum noise with (bin{t)bl^{t')) = 5{t — t') at zero temperature. The output field, describing photons escaping 
the cavity, is given by 5out = ^in + a and obtained by solving Eq. ( |T3] l. 

Squeezing properties of the cavity are revealed by homodyne detection of the output field ||5] |6l, mixing it with a local 
oscillator with the cavity frequency ujq and phase 9. The power spectrum iS'd(w) lH] of the homodyne detector signal /^(f) oc 
+ h.c. exhibits a deep at zero frequency and 9 = 7r/2 + '■p/2, where the squeezing effect is most visible, with 

= 1 +{A.Yf - 1) < 1. (14) 

bo K 

assuming weak external damping kq <C k. denotes the vacuum level, measured in the absence of ac excitation. A value 
smaller than one in this equation indicates a squeezed output field. Interestingly, output squeezing is weak in this limit due to 
pollution by the input vacuum noise. As discussed in appendix]^ the output field gets more squeezed as kq increases, the best 
squeezing being obtained for equal external and electronic dampings. This is however at the price of weaker squeezing in the 
cavity. Perfect squeezing could even be reached in the output field, limiting in this case the cavity state to half-squeezing. 


V. CONCLUSION AND OUTLOOK 


We studied the squeezing generated in a resonant cavity by coupling it to a mesoscopic conductor under parametric excitation. 
We showed that the quality of squeezing can be improved by enhancing nonlinearities in the conductor and by concentrating 
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the voltage excitation pulses at instants where the squeezed quadrature amplitude reaches its maximum. In optimal conditions, 
perfect squeezing can even be achieved. We remark that nonlinearities could also be enhanced in a tunnel junction by increasing 
the coupling to the cavity and the associated dynamical Coulomb blockade 1601 . More generally, our results can be easily 
extended to other quantum conductors for which the photo-assisted noise spectra are known or can be computed. They also 
suggest the possibility to engineer squeezed light for quantum information using electronic quantum conductors. 

Acknowledgments. We thank C. Altimiras, B. Huard, P. Joyez, F. Mallet, F. Portier, P. Simon and M. Trif for fruitful discus¬ 
sions. U.C.M. acknowledges the support from CNPq-Brazil (Project No. 229659/2013-6). 


Appendix A: Leviton excitation 


A Leviton is a pulse shape designed to minimize shot noise in driving electric current 1^ . In its periodized form, the 
time-dependent voltage is given by a sum of quantized Lorentzian pulses 


Vac{t) 


-2Hu;o + -y 

P 


l^Z 


1 r 

TT + {t — Itt/oJo)^ ’ 


(Al) 


of common width t. The ac phase is obtained by time-integration, (padt) = {e/h) J* dt'Vac{t'). It is conveniently written using 
the cyclic variable z = namely after a few algebraic manipulation 1551 


e 




1 z — a 
z 1 — az 


(A2) 


where 0 < a = < 1 is related to the pulses width. The conventional Leviton pulse, shaped to minimize zero-frequency 

shot noise, has the dc level eV = 2hujQ, cancelling the first term in Eq. ( |A1| |. In the main text, the dc bias chosen to optimize 
squeezing is eV = huJo. It can be understood by noting that the goal here is to minimize finite frequency noise correlators at 
the frequency ujq. Expanding Eq. ( |A2| l in powers of z, we obtain c_i = —a, c„>o = a"(l — a^): there is a strong imbalance 
between absorption and emission of photons, a = 1 is no longer an ac modulation but corresponds to a shift of the dc voltage 
by —2{h/e)ujo- 


Appendix B: Keldysh formulation and Heisenberg-Langevin equation 


The out-of-equilibrium physics of our system is conveniently described within the Keldysh path-integral formalism ED, 
enabling a systematic integration of the electronic degrees of freedom and yielding an effective action for the photons ll62l . The 
action obtained, Eqs. ( |B3| l and ( |B4| l, can be shown to be equivalent to the Heisenberg-Langevin equation of motion Eq. (|^ used 
in the main text, in the limit of small k/ojq <C 1. Eor a direct derivation of the Heisenberg-Langevin equation in the spirit of 
input-output theory, see Supplementary Note 2. 

The partition function (hereafter h = 1) 


Z = 



V[c, 


(Bl) 


involves an integration over complex-valued fields, a,a* for photons and c,c* for electrons. Sg is the action for the isolated 
conductor, S'cav the photonic part corresponding to the Hamiltonian iTcav = htoodPa and the electron-photon coupling is to first 
order 


5ep = 



(B2) 


where rj denotes the Keldysh time branch. I is the quantum conductor current written with complex-valued fields. 

The electronic part can be rigorously integrated using the cumulant expansion (e®'^'p)e = where we 

use the notations SSgp = S'ep — (^ep) and f I?[c, A = (A}g. We assume (S'ep)e = 0, a finite (S'ep)e can be absorbed by 

a small shift c>c p in the cavity fields. To summarize, we square Sgp, take its quantum average restricted to electronic degrees 
of freedom, and thus obtain a self-energy kernel for the photons involving current noise correlators. Eor clarity, we switch to 
classicaFquantum variables, Ud/q — (a+ ± a-)!s/2 and write the action in frequency space in order to take advantage of the 
current noise correlators given in Eq. Summing S'cav (i/2)(Sgp)e, we find the effective action 


Seff 





-Srr ) 



+ S< 


(B3) 
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where we used the notation and aci{t) = J^aci{uj)e The retarded photon Green’s function = 

a; — Wq + *k/ 2 has a pole shifted by half the damping rate k (see main text). The quantum-quantum self-energy part is = 
—iA^[5'o(a;o) -I- S'o(— wq)]. The effective action 5'eff includes also an anomalous term, responsible for state squeezing, 

Sa = -i\^ [ (a*(w)a*(2a;o - a;)S'i(a;o)-f c.c.) . (B4) 

J to 

Note that the real part of the self-energy induces in general a cavity pull which has been absorbed into a redefinition of loq. 
Computing this frequency shift consistently requires the second order term in the expansion of Ht in powers of g/oJo, namely 

2\2 r+ao 

=- 7^ dt'^f]VT,r,{t)[a*{t) - a^it)] , (B5) 


where Vr =T + T^ ■ The cavity frequency shift then vanishes for a tunnel junction. 

The derivation of Eqs. ( |B3| l and ( |B4| ) relies on the rotating-wave approximation, valid for k ujq, where the fields 
take significant values only for uj ~ ujq. In principle, the anomalous part of the action Sa also contains 
terms with aci{oj)aq{2uJo — w), corresponding to the effect of non-stationary noise terms on the damping. Those terms are found 
to be proportional to S'i(a;) — S'i(2a;o — w) and thus vanish for ui ~ wq, with the small parameter kJcoq. 

Finally, we discuss the connexion between the quadratic action of Eq. ( |B3| l and the Heisenberg-Langevin equation (|^. Quite 
generally, it is known that current fluctuations in a tunnel junction, or a quantum conductor, are not gaussian. Nevertheless, 
computing the non-gaussian current contributions to the statistics of photons, one finds that they are small compared to the 
dominant Wick-like contractions among the current operators. This is true in the limit of weak damping, n ujq, where cavity 
correlation functions only involve current operators I at frequencies iwQ. For example, the fourth-order cavity held correlator 
gives, for \ti\ 1/wo, 


(t)A(t)A(t)A(t)\ 


% 


Hi 



1 + 0 


ina,x{\eV/h- wo|,eKc/S) 


(B6) 


implying photon gaussian statistics except for the specific case of eV close to Huiq with no ac excitation. A related discussion 
can be found in Refs. El [32. The reason is that an electron-hole excitation with energy Hujq created by a current operator 
luiQ must be destroyed by another single current operator, the phase-space for alternative processes - where electron and hole 
are annihilated by two distinct current operators - being negligible for weak damping k. This argument pertains to higher-order 
correlation functions such that, for the purpose of photon statistics, it is legitimate to keep only the gaussian part of electronic 
current fluctuations. The resulting cavity held statistics are obviously gaussian. 

The above cumulant expansion can be rigorously stopped after the second order and the gaussian action in Eq. •HD becomes 
exact as long as k is negligible with respect to loq. Computing second-order cavity correlations functions, with different ordering 
of a and a^, we And coinciding results for the gaussian action and the Heisenberg-Langevin evolution. This completes the proof 
of the equivalence of the two formulations. This comparison differs from the standard derivation of a classical Langevin equation 
using the Keldysh action li62 . in which case information about operator ordering is lost. 


Appendix C: Optimized squeezing for a tunnel junction 


In this section, we set = 1 for simplicity. We focus on the zero-temperature case, relevant to maximize the cavity state 
compression. In this case. 


So [cV + (2n + Ijwo] — — sE -I- (2n + l)wo 

Hx L 


and the denominator in Eq. ([^ of the main text simplifies to 

(|c„p - |c„+ip) (eV + (2n -f Ijwo) = 2wo, 


(Cl) 


(C2) 


n^Z 


regardless of V and the c„ coefficients. At the optimal dc voltage eV = huiQ, the squeezed variance takes the simple form 

AXf = Y,\Cn±Cn-l\^\n\, (C3) 


n^Z 
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which we still need to minimize with respect to the distribution of Fourier coefficients c„, 

/ ‘^0 

We now prove that the pulse shape of Eq. (|^ in the main text, corresponding to the piecewise linear phase (^ac,opt(t) = 7r/2—wof, 
extremizes the variance AXf. We first differentiate Eq. ( |C4| l to obtain 

dlpacit) TT 

which gives —when evaluated at (/<ac,opt(i)- Using this result, we can proceed with the derivative of AXf 
with respect to an arbitrary form of (f) 



dAXf 

d'pe.cit) 


2^ |n|Re [« - - d^^^Cn-i)] ■ 


We evaluate this derivative with ^i>ac,opt(i) and its coefficients 

_ 1 1 
— I 7777) 

7rn + 1/2 

and obtain 


dAX^ 

9(/)ac(t) 


2wo 


Re 


^(n + l/2)(n-l/2) 


I ^—i(2n+l)uJot _ ^—i{2n—l)uJot' 


= 0 , 


(C6) 


(C7) 


(C8) 


which completes the proof. 

Inserting the coefficients Eq. (|C7|l into the quadrature variance Eq. (|C3|l, we find 


AXl 


1/2)2 


(C9) 


We also checked numerically that (j>a.c,opt{t) reaches the global minimum of AXj. 


Appendix D: Landauer-Biittiker calculation of the noise 

The noise properties of the quantum dot are derived using the scattering, or Landauer-Biittiker, formalism 0. The current 
operator is expanded over the basis of one-particle scattering states originating from both leads. The general case is reviewed 
in the Supplementary Note 4, we focus here on the asymmetric case where the probability for single-electron transmission at 
resonance AVjj/T <C 1 is small, and expressions simplify. Omitting spin, the current operator has the form 

^ zey^W f f deide2clj{ei)g{e2)cLie2) 

^ J J (Dl) 

y^d^[{ei-£2 + V)t/h+<l>^c{t)] _|_ ^ ^ 

with the Breit-Wigner resonant function g{e) = {e — £d + ifiT /2)“^. The operator cu/l{£) annihilates an electron in a scattering 
state of energy e incoming from the upper/lower lead. The normalisation is fixed by the average 

(c|j(£)ca'(e')) = (5a,a'/(e)<5(£ - s') (D2) 

with the Eermi function /(e) = (1 -f Due to the small capacitance at the upper dot-lead contact, the dc and ac bias 

voltages are applied essentially across this tunnel contact, the voltage potentials on both the quantum dot and the lower lead are 
fixed to the ground. Apart from the Breit-Wigner function, the rest of the calculation is similar to the case of a tunnel junction. 
The two-current correlators have the form of Eq. 0 and Eq. 0, where the equilibrium noise terms are given Eq. ( |TT] i in the 
main text. 
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Appendix E: Finite external damping 


We briefly discuss the case of a bare cavity damping kq comparable to the electronic damping k, but still much smaller than 
the resonator frequency wg- The complete Heisenberg-Langevin equation (El is solved by considering both the input held and 
electronic current fluctuations. One obtains for the two cavity held quadratures 


AX?/2(«o) = 1 + - l) , 

where 


AX2/2(0) = 


<S'o(u;o) + So{—ujo) 2Re [5'i(a;o)] 
5'o(<j-'o) — So{—ujo) 


(El) 


(E2) 


denote their variances in the absence of intrinsic damping kq, also given by Eq. (§ and discussed in length in the main text. 
Additionally, one finds for the output held squeezing, characterized by the power spectrum. 


Sd{oj = 0 ) 

5?, 


AkqK 
{k + ko)2 


(AA2(0) - 1) . 


(E3) 


Whereas a vanishing AAj(O) clearly optimizes squeezing in both the cavity and output fields, there is no such choice for kq. 
Increasing kq from zero improves squeezing in the output held but degrades cavity squeezing. Perfect squeezing in the output 
held is reached for k = kq, with vanishing AX^(O), in which case the cavity held is only half-squeezed. 
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CAVITY SQUEEZING BY A QUANTUM CONDUCTOR (SUPPLEMENTARY MATERIAL) 

Cited equations not preceded by S— refer to the main text. Unless indicated otherwise, we seth=l throughout the Supple¬ 
mentary material. 


S-I: Current-current correlation function 


In this section, we derive the current-current correlation function of an ac driven tunnel junction in the absence of cavity under 
the general form of Eq. (|^ in the main text. We use the definition 


/ OO 

-OO 

The current operator through the junction is given by 

(7kqe*[^*+'^“^*^'c|jk(t)cLq(f) - h-c.^ , 

kq 

corresponding to the operator 

^ c{jjj.CLq 
kq 


(s-1) 


(S-2) 


(S-3) 


in Eq. ([^ of the main text. Here cj^.^ (cka) is the electron creation (annihilation) operator in lead a = {U, L}, where U and L 
stand for upper and lower leads in a top-down geometry. In the absence of interaction and cavity photons, the Heisenberg time 
evolution of fermion operators is simply CLq(f) = Using the expansion in Eourier coefficients 


il)>ac(t) 


= E< 

riGZ 


2inuiot V 


(S-4) 


the two-time correlation function takes the form 


(/(fi)/(f2)) = l7kqp('c„i<^/(ek)/(eq) (S-5) 

kq ^ 

ni ,n2 





where we incorporated the phase terms in the Eourier coefficients Cne^^‘^ —>■ c„ and defined cukq = £k — £q- /(e) = (1 + ^ 

and f{e) = 1 — /(e) are the electron and hole Eermi-Dirac distribution functions. Substituting Eq. ( |S3] l in Eq. (|S3 and 
integrating over U and t 2 , we obtain 


(/(c.i)/a;2)) 



ni,n2 


dsf^dsq 


CmC*„J{ek)fiSq) 


X S{uJkq -b wi -I- eU -I- 2niUJo)6{uJi + uj 2 — 2(^2 — ni)a;o) 
+ Cn^Cn^f{Sk)f{Sq)S{Wkq+UJi - eV - 2niUJo) 

X S{oJi +UJ2 — 2 {ni — n2)uJo) . 


(S-6) 


In deriving this expression, we used the standard assumptions of weak momentum dependence of 7 kq and of the density of state 
in the leads, taken to be constant and denoted i/q. We introduced the tunnel resistance Rt = l/(27r|7pr'g). Before pursuing our 
derivation of current-current correlation function, we define the unsymmetrized equilibrium Johnson-Nyquist noise for a tunnel 
junction, in absence of both dc and ac bias, as 


Sioj) 


(/(f)/(0))e*‘^‘. 


(S-7) 


— OO 
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Using the expression Eq. (|S-2|i of the current operator, we readily hnd 


S{u}) = J dekdeqf{ek)f{eq)6{uji^q + u}) 

1 2a; 


Rt (1 - 

Using the above dehnition, we are able to rewrite the current-current correlator [Eq. S] as 

{i{uji)iuj 2 )) = ^ 5'„(a;i)27ri5(a;i + uj 2 - 2na;o), 


(S-8) 


(S-9) 


where the different noise terms are dehned as 

-S'n(wi) = ^ X! [cn'C*/+„5(a;i -I- eU -I- 2n'a;o) 

n-i 

+ C*^,Cn'-nS(LOi- eV -2n'LOQ) . 

We thus obtain Eqs. Q and Q of the main text once the phase terms c„ —>■ have been reinstated. 


(S-10) 


S-II: Heisenberg-Langevin equation 

We derive in this section the Heisenberg-Langevin equation, Eq. (|^ in the main text, following the steps of standard input- 
output theory m. A similar analysis can be found in Ref. 12. The cavity is coupled to a quantum conductor and also to external 
electromagnetic modes. Both play the role of an environment for the cavity, absorbing and emitting photons via quantum vacuum 
noise terms. The Hamiltonian of the complete system is 

H = LOQb)a + ^ eakcjjkCak + “ d)i 

ak 

-aS|j). (S-1) 

n n 

The first term describes the isolated cavity, the second one the scattering states c^k El incoming from both leads a = U / L, 
and which includes scattering by the quantum conductor. The third term is the coupling between the cavity single-mode and 
the electric current through the conductor. Eor a tunnel junction, the scattering states Cak coincide with lead electrons to lowest 
order in the tunneling. The current operator I in the Hamiltonian is then 

i = i^ (^7kqe*[^*+‘^“(‘^lc|jkCLq - h.c.) . (S-2) 

kq 

It is in the Schrodinger picture where the only time dependence is through the voltage. In the rest of this section, we will focus 
our analysis on the tunnel junction for simplicity. 

The last two terms in Eq. ( |S-1| ) stand for the external photonic modes whose coupling to the cavity is assumed to be weak. We 
now switch to the Heisenberg picture where d{t) = e*^*ae“*^* and compute the equation of motion (EOM) for the cavity field 
from d{t) = i\H, d] (t). The external modes can be treated perturbatively, as detailed in Ref. [[T], reducing the EOM to 

d{t) = -iu!od{t) - + Xlnit), (S-3) 

where kq denotes the cavity damping rate due to these modes. 6in(f) is an input held which feeds the cavity with vacuum noise 
at zero temperature. The cavity held is also coupled to the output held (from which cavity properties can be extracted) through 
the boundary condition 


&out(f) = hn{t) + y/i^d{t). (S-4) 

In{t) in Eq. (F2 denotes the current in the Heisenberg picture. It is evolved with the complete Hamiltonian H including, via 
coupling to the cavity, the cavity single-mode and its electromagnetic environment. We will write the time evolution of dait) 
perturbatively in the cavity-conductor coupling A and insert the result into the cavity EOM. 
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Let us start by introducing the Heisenberg-picture operators 

/kq(0 = c[q(0cuk(i) (S-5) 

which appear in the current Init)- Taking the derivative with respect to time, one obtains 

fkq{t) = -*(£Uk - eLq)/kq(i) ” A[Vi„t (t) , /kq (t)] (S-6) 

in which ~ Assuming a weak cavity-conductor coupling, we solve Eq. ( |S-6| l perturbatively. The 

lowest order in A is readily integrated to give 


/kq(0 = /kqe-*‘^-‘o‘, (S-7) 

where we define Wkq = euk ~ SLq- At this order, the time evolution does not include the cavity mode and the current operator 
is given by Eq. ( |S-2[ ) and is denoted I{t). Replacing /j^q(f) in the last term of Eq. ( |S-6[ ) by the zeroth order solution Eq. ( |S-7| i, 
and the current by I{t) in Vint{f), we obtain the first order correction 


/kq(0 = - A[Vint(f), /kq]e 

This differential equation can be integrated with the result 


/kq(i) = /kqe-^^-'o* -X f 

J —OO 


(S-8) 


(S-9) 


We consider the current and obtain 


/H(f) = _ h.c.) 

kq 

= Ht) - A /" [a^(r) - a(r)][/(T), 7 (f)]dr. 

J —OO 

Substituting Eq. ( |S-10| l in Eq. ( |S-3| l we obtain 

a{t) = -iiOQaiyt) - ya(f) - y^&in(f) + A/(f) 

-A^/ [a^r) - a{T)][i{T)J{t)]dT. 

J —OO 


(S-10) 


(S-11) 


The last term in this expression is second order in A. We thus average it with respect to the quantum conductor Hamiltonian 
isolated from the cavity, and use the expansion Eq. ( |S-91 l to rewrite 

([J(r),7(f)]) = ^ [ 

"ez"' (S-12) 

^ ^^iu}(t—T) — 2inujQt —r) —2maJoT^ 

Our approach assumes small kq and A such that the cavity field is essentially oscillating at the frequency wq. After the change 
of variable 


a{t) (S-13) 

d{t) becomes a slow field which evolves over time scales much larger than the decay of [7(r), 7(f)]. We therefore replace (t) 
by (f) in the last term of Eq. ( |S-11| ) which becomes 

m = -yaW - + A7(f)e“«‘ 

+ A^C'_|_(f)a’^(f) — A^C_(f)o(f), 


(S-14) 






16 


with 


C±(i)=r dre*‘"“(‘±")([/(t),/(T)]). (S-15) 

J —OO 

Computing C± with Eq. ( |S-12| l, we find an imaginary principal part term, describing frequency pull of the cavity, which we 
absorb into the cavity resonant frequency. We apply the rotating-wave approximation to the real part discarding all fast oscillating 
terms (beating with frequencies multiple of wq) and obtain 

ReC+ ~ i[^i(a;o) - 5i(wo)] = 0 (S-16a) 

ReC_ ~ i[S'o(wo) - S'o(-wo)]- (S-16b) 

We use these results in Eq. ( |S-14| l and go back to the original time frame for the cavity field a{t) —?> to find the 

Heisenberg-Langevin equation 

a{t) = -iuJod{t) - ^ - v^^in(f) + A/(f), (S-17) 

advertised as Eq. ( [T3] l in the main text, where we introduced the cavity damping rate k = 2A^Re C-. In the absence of external 
electromagnetic modes kq = 0, it reduces to Eq. (|^ in the main text where the cavity only interacts with the quantum conductor. 

The solution to this equation can be obtained in time or frequency space. Combined with the boundary condition Eq. ( |S-4| i, it 
leads to an expression for the output field in terms of the input field and the current in the quantum conductor 




[w — Wo + i{K — Ko)/ 2] 4n(w) + iXy/K^I{uj) 
uj — ujQ + i(K + ko )/2 


(S-18) 


S-III: Finite-frequency Noise of a Quantum dot under ac excitation 


We consider a single-level quantum dot attached to two leads, denoted upper and lower leads, with tunnel couplings tjj and 
t]^. Ignoring weak Coulomb interaction on the dot and assuming spinless electrons, the Hamiltonian takes the form 



(S-1) 


with the constant density of states vq for the two leads. Ed is the dot single-level energy, annihilates an electron with energy 
e in lead a = U/L and satisfies the anticommutation relation I?) 

{cl(e)cQ'(e')} = Sa,a'S{e - e'). (S-2) 

The transport of electrons in this system is described by using scattering states. The scattering matrix is given by H 


fSuuis) Sul{£)\ 
\Slu{£) Sll{£)J 


-l + ig{s) 


( 


rr, ) 


(S-3) 


with the function g{e) = {e — Ed + /2) We have introduced the dot escape rates towards the two leads Tq, = Tri/olfaP, 

and the total escape rate r = rL + ru. We write the current operator in the presence of an ac driving of the two leads following 
Ref Q, where the scattering problem in the dot is disentangled from the potential modulation and the photo-excitation in the 
leads ©. Operators c^(e) are introduced, which describe scattering modes just after entering or leaving the dot. They are 
assumed to be adiabatically connected to the lead scattering modes through 


Ca(e) ='^Co{£- ‘2rnuJo)c^ 

m^Z 


(S-4) 


where we introduced the Fourier coefficients for lead a 


giej* dt' VAC,c{t') 


E ol 2 imuiot 


(S-5) 
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Eq. ( |S-4| i describes the fact that an electron with energy e arriving at the dot-lead contact is composed by a coherent sum of lead 
electrons with energies e — 2mujQ which have absorbed or emitted m photons from the ac field with probability amplitudes c^. 
Writting the current in terms of operators c^(e), as if no ac modulation were present, and substituting with Eq. ( |S-4| l, one obtains 
the general expression for the current operator in lead a 


Lit) = lf f dside, Y. 

n\,n 2 7:7^ 

xcj^(ei - 2niuiQ)A^^>{a,ei,e2)c^'{e2 - 2n2Wo), 

where we use the scattering amplitudes 

(g;, , £ 2 ) - Syy'Scfy' (S 1 ) (^2 ) 

in terms of Eq. Einally, the lead scattering operators have thermal distributions 

(cJ,(e)cQ/(e')) = da,a'fai£)Sis - e'), 


(S-6) 


(S-7) 


(S-8) 


where the dc voltage Va sets the chemical potential in lead a, namely fa{e) = (1 -f vA/^bT'^ i jsjQjg always 

move the dc voltages Va into the time dependence of Eq. (|S-6|), instead of having them in the lead thermal distributions. 


1. Asymmetric quantum dot 

Expressions simplify in the limit of an asymmetric dot, Ec/ E^, more strongly coupled to the lower lead than to the upper 
lead. Computing the following scattering probabilities 

|4lLL(t/,ei,£2)P = r^r2 \g{e,)\^\g{e2)? 

\Auu{UtSi,£2)\’^ = \All{U,£i,£2)\^ 

-f 4^sin2[5(ei) -(5(e2)], 

where we defined the scattering phase 

5{e) = arctan[E/(2(£d — £))], 

we see that they both scale as (Ej//E)^. In contrast to that, for E;/ <C E^, we find 

Aul{U,£ 1 , 62 ) — ig{£ 2 )\/^u^, 

Alu{U,£i,£2) — —ig*{£l)V^U^, 

and their square modulus scale as Efj/E. We therefore neglect the LL and UU terms in the current. In addition, we assume that 
the voltage drop occurs at the upper dot-lead contact so that = <5^.0 we note = Cm- The dc voltage is V = Vjj at 
the upper lead and Vl = 0 in the lower lead. The potential of the dot is also to the ground. With these different notations and 
approximations, the current lu (f) is given by Eq. in the main text (see Methods) where the dc voltage is included in the 
operator time-dependence. 


(S-9a) 

(S-9b) 


(S-lOa) 

(S-lOb) 


2. Fully symmetric quantum dot 

We briefly discuss the case of a quantum dot symmetrically coupled to the leads E[/ = Ej^ and to the cavity. The first 
order cavity-dot then takes the form iX{Iu — lL){dX — d) involving the upper and lower currents with equal weights. The 
general approach devised for the tunnel junction also applies in this case, with the current I being replaced by lu — Il in the 
Heisenberg-Langevin equation (Eq. in the main text). Setting 

I = iu- IL 


(S-11) 
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and computing the correlators {I{uji)I{uj 2 )), we obtain 


5'o(±wo) 


2 'Kh? 


de V 


^71 ^ 72 *^ 71*^72 
^ni^n2 ^ns '^ 713 + 722 —ni 


niTj.2'n.3 

7l'>'2 

^ ^7172 ^(^3 ~ ^ i Wo + 2(713 ~ '^l)wo] 

X ^7172 (e, e ± wo)/^! (e - 2 nia;o) 

X - 2n2Wo), 


(S-12) 


1 „72*„72 „7l* 

•''n2 ni+"3-’22 + l 


with /(e) = 1 — /(e) for the hole distribution, and 

^ 711712^3 

X 24 ^ 2 ^,! [e + Wo + 2(773 - 77 ' 2 )wo, e + 2(773 - 772 + l)wo] 

X A^j^2 (^> £ + (e - 277 iwo) 

X J ^ A ^ + ^O - 2772Wo). 

We also introduced the scattering amplitudes 

^7172(^11^2) = 24-^172 (C^J £lj £2) — 24 ,yj.y 2 (1^1 £ll £2). 


(S-13) 


(S-14) 


The knowledge of these correlators is sufficient to compute the quadrature variances of the cavity field, namely 


AXf/2 = 


‘S'o(wo) + <S'o(~<x'o) T 2Re [S'i(wo)] 
< 50 (^ 0 ) — S'o(—Wo) 


(S-15) 


We now turn to the numerical minimization of the variance AXf. The voltage is chosen to be zero on the dot so that the voltage 
potentials in the leads are opposite Vu{t) = —VAt) by symmetry (the voltage drop is the same at each lead-dot contact). The 
applied voltage has a constant dc part and a sinusoidal part, Vu{t) = V + Vi cos(2wof), such that the c„ coefficients are given 
by Bessel functions 


At a general level, squeezing is optimized by setting T = 0 and by having a sharp dot resonance T/wq <C 1. In this fully 
symmetric situation however, the condition eV = Hujo is no longer met to optimize squeezing. Focusing on the asymptotic limit 
r/wo 0, we find the best squeezing 


AX? = 0.654 for = -0.58, 

2/kJo 

ffd = 0 and eV = O.Gwq. Squeezing is in fact not very sensitive to V and Sd for F ^ wq as long as |I4| > 


(S-17) 
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